Spectral features due to inter-Landau-level transitions in the Raman spectrum of 

bilayer graphene 
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We investigate the contribution of the low-energy electronic excitations towards the Raman spec- 
trum of bilayer graphene for the incoming photon energy > leV. Starting with the four-band 
tight-binding model, we derive an effective scattering amplitude that can be incorporated into the 
commonly used two-band approximation. Due to the influence of the high-energy bands, this ef- 
fective scattering amplitude is different from the contact interaction amplitude obtained within the 
two-band model alone. We then calculate the spectral density of the inelastic light scattering ac- 
companied by the excitation of electron-hole pairs in bilayer graphene. In the absence of a magnetic 
field, due to the parabolic dispersion of the low-energy bands in a bilayer crystal, this contribution 
is constant and in doped structures has a threshold at twice the Fermi energy. In an external mag- 
netic field, the dominant Raman-active modes are the n~ — >■ n"*" inter-Landau-level transitions with 
crossed polarisation of in/out photons. We estimate the quantum efficiency of a single n'^ 
transition in the magnetic field of lOT as ~ 

PACS numbers: 73.22.Pr, 71.70.Di, 78.67.Wj, 81.05.ue 



I. INTRODUCTION 

Bilayer grapheneii^ is a representative of the re- 
cently discovered family of new carbon allotropes.^ It 
attracted attention by the observation of an unusual se- 
quencing of plateaus in the quantum Hall effect^'^ and 
the possibility to modify its spectrum by opening a 
small gap and induce an insulating state with an exter- 
nal electric fieldi^"— The electronic properties of bilayer 
graphene have been characterised using angle-resolved 
photoemission spectroscopj*^ and optical absorption in 
the visible^'-- and infrared^^JA spectral range. Bilayer 
graphene has also been studied using inelastic scattering 
of light, subject to the detection of Raman-active lattice 
vibrations in this two-dimensional crystalJ^^— However, 
no theoretical or experimental study has yet been per- 
formed, on the contribution of electronic excitations to- 
wards Raman spectra of bilayer graphene. This paper 
presents a theory of such a contribution. 

The electronic Raman spectroscopy can provide infor- 
mation about various single particle and collective elec- 
tron excitations in the system studied. In semiconduc- 
tors, it has been, for example, employed to investigate 
donor and acceptor states, plasmons and spin-density 
fluctuations involving electron spin-flip due to the spin- 
orbit interactionj2^i2^ The inelastic scattering of photons 
on electrons in semiconductor placed in an external mag- 
netic field was first discussed in Ref. [2^ where it was 
pointed out that nonparabolicity of the electronic bands 
is crucial for the electron-photon interaction matrix ele- 
ments not to vanish. The features corresponding to the 
electronic contribution to the Raman scattering in an 
external magnetic field were observed in many semicon- 
ductors, for example, InSb (Ref. ^2f) and GaAs.^ 

Recently, the Raman spectroscopy of electronic ex- 
citations in monolayer graphene has been investigated 
theoreticallyi^ It has been shown that at high mag- 



netic fields the inelastic light scattering accompanied by 
the excitation of the electronic mode with the highest 
quantum efficiency involves the generation of inter-band 
electron-hole pairs. At high (quantizing) magnetic fields 
this leads to the electron excitations from the Landau 
level (LL) n~ at the energy —\/2nhv/XB in the valence 



band to the Landau level n'^ at the energy \/2nhv/ \b 



in the conduction band with energies aj„ = 2\/2nhv/ Xb 
and crossed polarisation of in/out photons, in contrast 
to the An — ±1 transitions between Landau levels 
which are dominant in the absorption of left and right- 
handed circularly polarised infrared photons,^'' Raman 
spectroscopy, therefore, provides data supplementary to 
that obtained in optical absorption. This fact could be 
of interest in particular for bilayer graphene placed in an 
external magnetic field. Recent measurement of the in- 
frared absorption in fields up to 20T showed^^ deviations 
from the tight-binding model for an ideal bilayer. Here, 
we offer a theory for an alternative experimental probe, 
with different selection rules for the inter-Landau-level 
excitations, which could bring some new insight into the 
properties of the bilayer graphene. 

We study the Raman spectroscopy of electronic exci- 
tations in bilayer graphene both with and without an 
external magnetic field using the tight-binding approach. 
First, we describe processes in which after the inelastic 
scattering of the optical photon, an electron-hole pair is 
created in the low-energy bands of the bilayer. Then, 
we evaluate the scattering amplitude corresponding to 
such a process within the framework of the four-band 
tight-binding model and extract an effective scaterring 
amplitude that is incorporated into the two-band low- 
energy Hamiltonian for bilayer graphenej^ This effective 
scattering amplitude contains the influence of the high- 
energy bands and is different from the scattering ampli- 
tude found using only the approximate low-energy Hamil- 
tonian. Next, we find the low-energy electronic contribu- 
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tion to the Raman spectra and obtain an analytic descrip- 
tion for the selection rules and intensity of such Raman 
scattering in a bilayer placed in a strong (quantizing) 
magnetic field. 

II. THEORY OF THE INELASTIC LIGHT 
SCATTERING IN BILAYER GRAPHENE 

Bilayer graphene consists of two coupled sheets of 
graphene with AB (Bernal) stacking characteristic of 
bulk graphite^ Fig. 1(a). The unit cell contains four in- 
equivalent atoms Al, Bl, A2 and B2 where letters A and 
B denote two sublattices in the same layer and 1/2 stands 
for the bottom/top layer. The Fermi level in graphene 
lies in the vicinities of the corners of the hexagonal Bril- 
louin zone (also called valleys) known as and K- 
(Fig. 1(b)). The conventional tight-binding Hamiltonian 
based on 7r-orbitals of carbon atoms (one per atom, four 
in the unit cell) and expanded in momentum around the 
valleys reads 

H = ( ^^3{<JxPx - CTyPy) £.VCrp\ 

Here, cr — {ax,<Jy) and ax,Oy,az are the Pauli matrices, 
electron momentum p is measured from the center of the 
valley, v ~ lO^m/s (Refs. UUIIilillil) is a parameter 
related to the nearest neighbour intralayer c ouplin g 70 ~ 
3gY n.^^^ 71 ^ 0.4eV (Refs. ilMUlillilllllll) is the 
direct interlayer coupling, is related to the weak direct 
Al <-> B2 interlayer hops [v^/v ~ 0.1 (Refs. Qliilll)] and 
^ = ± is the valley index. The basis is constructed using 
components corresponding to atomic sites Al, B2, A2, Bl 
in the valley K+ and B2,A1,B1,A2 in K^. One can 
also take into account terms quadratic in the electron 
momentum p: 

\ ^x{pl-pl)-'^(JyPxPy 
2 

where /x = — However, the influence of the 5H term 
on the results of the Raman spectra analysis is negligibly 
small, as shown in the Appendix. 

The part of the resulting electronic dispersion relevant 
for the Raman scattering of photons with energies Q, < 
2.5eV, is illustrated in Fig. 1(c) for the valley K+. Two 
bands, later referred to as low-energy ones, touch each 
other at the neutrality point - the position of the Fermi 
energy in the neutral structure. Two other, referred to as 
high-energy bands, are split by the interlayer coupling, 
71, from the neutrality point. The U3 parameter leads 
to the trigonal warping of the electronic dispersion. Its 
influence is most important for very low energies, e < 
5meV. 

To describe the process of inelastic scattering of light 
on electrons, we consider an experimental setup in which 
incoming laser light of energy f2 3> 71 , in-plane momen- 
tum q (out-of-plane component of momentum equal to 




FIG. 1: (a) Schematic of the bilayer graphene crystal lattice, 
(b) The Brillouin zone of bilayer graphene with two inequiv- 
alent valleys denoted as K+ and K- . (c) The band structure 
of bilayer graphene in the vicinity of the K+ point along the 
Px axis. Also shown are some of the two-step processes lead- 
ing to the creation of an electron-hole pair in the low-energy 
bands accompanied by the absorption of a photon followed by 
emission. Gray solid (dashed) lines indicate the first (second) 
step of the process. The black (white) circle denotes the hole 
(electron) in the final electron-hole pair, while the hatched 
circle represents the intermediate virtual state. Note that for 
any intermediate state \u) with energy e^, 
Diagrammatic representation of the scattering amplitudes TZ, 
discussed in Sec. II, and SR. (discussed in Appendix). 



Qz = y^Q'^/c'^ — q^) and polarisation 1 is shined onto to 
the sample. Scattered photon has polarisation 1, in-plane 
momentum q and energy = fl — uj, where uj is the 
Raman shift. We also assume the temperature T to be 
smaller than the Raman shift, fc^T < w (fc^ is the Boltz- 
mann's constant). In our case, the inelastic light scatter- 
ing may occur via a one-step process (so called contact 
interaction) or a two-step process involving an interme- 
diate state. The two-step process, such as shown in Fig. 
1(c), involves: the absorption (or emission) of a photon 
with energy il (fl) transferring an electron with momen- 
tum p from an occupied state in the valence band into 
a virtual intermediate state (energy is not conserved at 
this stage) , followed by another electron emission (or ab- 
sorption) of the second photon with energy f2 (il). The 
one-step process is the usual inelastic scattering of an 
incoming photon on an electron with transfer of energy 
to the latter. As a result of both one and two-step pro- 
cesses, an electron-hole pair in the low-energy bands is 
created with the electron and the hole having almost the 
same momentum (p + q — q and p, respectively), since 
q, q ^ p and the momentum transfer from light is neg- 
ligible {v/c ~ 3 • 10"'^). Therefore, p -|- q — q « p and 
due to the approximately electron-hole symmetric band 
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structure in the vicinity of Brillouin zone corners, the 
electron initial and final energies and are related, 

To include the interaction of the electrons with 
photons, we construct the canonical momentum [p — 
e(A(r, t') + A(r, t"))] , where A(r, t') and A(r, t") are the 
vector potentials of the incoming and outgoing light, re- 
spectively, 



h.c 



(2) 



and bq,5j,i is an annihilation operator for a photon with 
in-plane momentum q, out-of-plane momentum compo- 
nent qz and polarisation 1. We expand the resulting 
Hamiltonian up to the second order in the vector po- 
tential and write down the interaction part, 

Ai„t = j - (A(r, t') + A(r, t")) + ^ ^ , (3) 



2 ^ dpidpj 



,dHo 
dp 



where j : 

We aim to calculate the spectral density .g(w) and the 
quantum efficiency (intensity) of the Raman scattering, 
/. The quantum efficiency describes the ratio of the flux 
of outgoing, inelastically scattered photons to the flux of 



the incoming photons, and is an integral, / = Jdujg{u}), of 
the spectral density g{uj) representing the probability for 
the incoming photon to scatter inelastically with energy 
= — cj, where uj is the Raman shift. 

The quantum efficiency expresses the total probabil- 
ity for single incoming photon to scatter inelastically in 
a proccess under consideration - that is, to scatter on 
an electron and excite an electron-hole pair in the low- 
energy bands. The probability for the incoming photon 
to scatter with the Raman shift w in a particular direction 
(defined by the momentum q of the scattered photon), 
is, in turn, characterised by the angle- resolved probabil- 
ity of scattering u'(q). Finally, the scattering probability 
w(q) that one photon is scattered with the excitation of 
an electron-hole (e-h) pair in the final state is related, as 
w oc |7?,p, to the scattering amplitude TZ of the Raman 
process. 

The amplitude TZ is the sum of the amplitudes corre- 
sponding to the one-step and two-step processes. Only 
terms quadratic in electron momentum p which appear 
in the addition SH to the Hamiltonian in ([T|), contribute 
to the contact interaction. We show in Appendix that 
this contribution is much smaller than the leading con- 
tribution from the two-step processes, thus, we neglect 
it in further considerations. To find TZ, illustrated us- 
ing Feynman diagrams shown in Fig. 1(d), we describe 
a two-step transition which involves an intermediate vir- 
tual state 1 1/) with energy e,y, as 



7^ = - 




ei(^/-^-)*'(j.r)e-^(^--^*')|i/)(zy|e^('i--"*")(j-l)e^(^--^-)*"di'di"- 

ei(^/-^-)*'(j.l)et('i-"--"*')|,.)(iy|e--^('i----^*'')(j-r)e*(^''-^-)*''di'dt''. 

I 



(4) 



The virtual state may belong to any of the four bands, 
since an electron is excited from a state with momentum 
p to a state with momentum p -I- q or p — q depending 
on the accompanying photon process. At this step of the 
calculation we still work with the four-band Hamiltonian 
([1]), to include the influence of the high-energy "split" 
bands. In Eq. (U), the flrst (second) term corresponds to 
processes in which the photon is absorbed (emitted) in 
the first step and emitted (absorbed) in the second step 
of the process and is given by the first (second) diagram 
in the expression for TZ in Fig. 1(d). Integration in the 
time-dependent perturbation theory in Eq. can be 
performed by changing variables to t = i' — i", which 
varies at the scale of w"^, oj = fi — Q., and t ~ {t' + t")/2, 
which varies at the scale of ^ = {il+n')/2 ^ lu. For 
incoming and outgoing photons, 3> 71, and we also 
study the low energy excitations in the final states with 



cj ^ 71 . This allows us to expand factors resulting 
from the integration over r in powers of (ciz/fi), keeping 
terms of the order of 1 and (71 /Q) [the latter appear when 
the virtual state is taken to be in the high-energy bands] 
and to perform summation over the intermediate virtual 
states of the process. Consequently, the amplitude TZ 
takes the form of a matrix 



TZ 



<Jz 



(ixr)^ 



M d 



eon^ \ "K az J\ Jz ' n 
jiay ^v{ayPx + (JxPy) 



(5(e/-ej-cj); 



^v{o-yPx + CTxPy) 







llCTx ^v{(TxPx - O-yPy) 



Below, we analyse the contribution of electronic modes 
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toward the low-energy part of Raman spectrum with the 
photon energy shift to < 71 /2, which is determined by the 
excitation of the electron-hole pairs in the low-energy (de- 
generate) bands with ^ 71. At such low energies, the 
band structure as well as Landau level structure can be 
described by the effective two-band Hamiltonian written 



Hcs = -— [{pI-pD Ox + . (5) 

To describe the excitation of the low-energy modes cor- 
responding to the transitions between low-energy band 
states described by iJcff', we take only the part of 7?. which 
acts in that two-dimensional Hilbert space, keep terms in 
the lowest relevant order in vp/71 ^ 1 and 71 /ri ^ 1, 
and write down an effective amplitude 7?.off, 



(Tyd^]y (6) 



We point out that the above matrix cannot be obtained 
within a theory constrained by the two-band approxi- 
mation, Eq. ([5]), from the very beginning. Seemingly, 
one may try to define a contact-interaction-like term due 
to the terms quadratic in the electron momentum p in 
Eq. (|5j), which carries a prefactor ^^^p^, which may 
suggest a greater magnitude of scattering than prefac- 

2^2 2 

tor ^^^2 above. However, the scattering amplitude ob- 
tainecf within this model can only be applied to photons 
with ri < 71, which is hardly relevant for Raman spec- 
troscopy since the latter is usually performed with laser 
beams using ~ 1.3 - 2.8eV^"^ 

The angle-resolved probability of the Raman scatter- 
ing, w{ci ~ 0), determined using Fermi's golden rule and 
with the help of Eq. ^ , is 



rfp|(/|7^cff|^)P X /, (1 - ff) d{e,+LO- ej) , 



CTT 



ft3 



where fi and // are filling factors of the initial and fi- 
nal electronic state, respectively, and the spin and valley 
degeneracies have already been taken into account. This 



givei 



,33 



Tie^/M^ r„ 7? „ , 



(7) 



_ 2 
1 X 1* 



1 + (1 X r) • (i X i* 



Above, the first term with polarization factor Sg de- 
scribes the contribution of photons scattered with the 
same circular polarization as the incoming beam. The 
second term, with polarization factor So, represents the 
scattered photons with circular polarization opposite to 
the incoming beam. 

In turn, the angle-integrated spectral density of Raman 



scattering g{oj) is 



(8) 



Here, the constant spectral density as a function of 
Lj reflects the parabolicity of the low-energy bands and 
thus, energy-independent density of states in the bilayer. 
This is different in monolayer graphene, where g{uj) oc 
w, reflecting the energy-dependent density of states of 
electron-hole pairs. The characteristic of monolayer 
graphene crossed polarisation of in/out photons is re- 
tained in the case of the bilayer system. Experimentally, 
constant spectral density g in undoped bilayer graphene 
is impossible to distinguish from a homogeneous back- 
ground. However, if the chemical potential /i is not at the 
neutrality point, then transitions with lo < 2fj, are essen- 
tially blocked. Although new processes, resulting in the 
creation of the intraband electron-hole pair excitations 
and very small w, are possible for /i 7^ 0, their contribu- 
tion carries additional prefactor w/c ~ ^ Explicit cal- 
culation performed for the monolayer graphene showed 
that the quantum efficiency of the intraband transitions 
was of the order of 10~^^^ In contrast, for chemical po- 
tential /I ~ 50meV (corresponding to additional carrier 
density no ^ 1.5 x lO^^cm"^), the lost quantum efficiency 
due to the blocked interband transitions is, according to 
Eq. dll), A/~ 10-12. 



III. INTER-LANDAU-LEVEL TRANSITIONS IN 
BILAYER GRAPHENE RAMAN 

The quantization of electron states into Landau levels 
gives the Raman spectrum due to the electronic excita- 
tions a pronounced structure which can be used to detect 
their contribution experimentally. We only consider here 
low-energy Landau levels, as at high energies the Landau 
level broadening due to, for example, electron-phonon in- 
teraction, will smear out the LL spectrum. In strong 
magnetic fields, low-energy Landau levels are sufficiently 
described'^'* by 



71 A| 

V'n 





y/n{n~ 1); 
, n = 0, 1; 



(9) 



n> 2: 



where — yhJeB is the magnetic length, n is the 
Landau level index and a = -f denotes the conduction 
and a = — the valence band. Also, V'n is the normalised 
n-th Landau level wavefunction. In a neutral bilayer, 
all LLs have additional fourfold degeneracy (two due to 
the electron spin and two due to the valley). Moreover, 
levels n = and n = 1 are degenerate at e = giving rise 
to an 8-fold degenerate LL. We can project our effective 
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transition amplitude 7?.eff onto the eigenstates \l/„c» to 
find the electronic Raman spectrum in the presence of a 
strong external magnetic field. This leads to the following 
selection rules for allowed electronic transitions from the 
initial level n~: 



Raman shift 



i) n 



ii) (ti =F 1) ^ ± 1)" 



(10) 



FIG. 2: (a) Schematic of allowed inter-LL transitions accom- 
panying the Raman scattering. The solid (dashed) line rep- 
resents the first dominant (weaker) transition 2~ <~ 2^ (pair 
2~ <— and 0-<— 2^). (b) The low-energy electronic contribu- 
tion to the Raman spectrum in bilayer graphene. The solid 
(dashed) line represents the spectrum in the presence (ab- 
sence) of an external magnetic field and chemical potential 
= (/^ 7^ 0). For the spectrum in a magnetic field, cor- 
responding inter-LL transitions have been attributed to each 
peak. 



Among those, i) is the dominant transition. These selec- 
tion rules, represented schematically in Fig. 2(a), show 
that using Raman spectroscopy, one can probe different 
electronic excitations than in optical spectroscopy, where 
the selection rules are An = ±l^i^ For a neutral bi- 
layer, the angle- integrated spectral density g{ijo) of Ra- 
man scattering in the magnetic field is equal to: 



g{uj) « 16^, 



hv 



4:TTeohc cj \Xb^ 



n>2 



2 \ 2 / fc \ 2 

e V \ hv 



^7(w-2e„+) + &/(o^) 



(11) 



^~e(K+i)+) + y^7(^--g(»+i)+ -£(«-!)+, 

n=l,2 ri>3 



Here, we use Lorentzian 7(2;) — 7r^^r/(x^ -f F^) with a 
width specified by F to model the broadening of Landau 
levels. The term 5g{ijj) describes the spectral density of 
the (n =p 1)^ — !■ (n ± 1) transitions, which is a correction 
to the dominant contribution due to the n~ — >■ tran- 
sitions given by the first term on the right hand side of 
Eq. (HH). 

An example of the low-energy electronic contribution 
to the Raman spectrum in the neutral bilayer in strong 
magnetic field is shown with a solid line in Fig. 2(b). The 
dominant features are peaks due to the n~ — >■ transi- 
tions with the first being the 2" — )■ 2+ transition. Note 
that within the LL indexing scheme applied here, indices 
and 1 are only used to denote one valley-degenerate 
level each (no a index is needed). Lifting the valley de- 
generacy by introducing charge asymmetry between lay- 
ers will not allow any additional n~ — ^ n"*" transition be- 
cause valley-split levels for n — 0,1 belong to different 
valleys and excitation between them would require a huge 
momentum transfer. The quantum efficiency of a single 
>n+ peak in Fig. 2(b) is approximately 



per incoming photon, which at the field B ^ lOT gives 
/„-_;.„+ ~ 10~^^ for J7 ~ leV photons, comparable to 
similar transitions in monolayer graphene*?^ 

A weaker feature in Fig. 2(b) is the first and the only 
visible (n=F 1)^ (n ± 1)^ peak due to both 2^ — )■ and 
— >■ 2+ transitions, positioned to the left of the 2~ — >-2+ 
peak. The quantum efficiencies of the (n±l)~— >-(n=Fl)+ 

transitions are smaller by the factor (jf-)^ in comparison 
to the n~ — transitions. This is different from the 
monolayer graphene case, where the corresponding ratio 
between quantum efficiencies of (n ± 1)^ — > (n =F 1)+ and 

transitions is (^) , much smaller than for the 
bilayer. The term 5g{uj) can be further emphasized by 
changing the energy of incoming photons fi. Shown in 
Figs. 3(a) and 3(b), is a comparison of the total spectral 
density g{uj) and contributions due to each mode sep- 
arately, for two different energies of incoming photons, 
= 2eV and = leV. The intensity scale is the same 
on both figures and in each case, the total spectral den- 
sity (?(w), the contributions due to the n+ and 
(n ± 1)~ — >■ (n =F l)"*" modes are shown in the solid, dot- 
dashed and dashed line, respectively. The dominant con- 
tribution, resulting from the Raman scattering accom- 
panied by the n~ — > electronic transitions, is pro- 
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■J 2 2J2 2-J6 a/ a, -Jl 2j2 2^6 aJfo, 

Raman shift Raman shift 



FIG. 3: Comparison of electronic contributions to the Raman 
spectra in neutral bilayer graphene for two different energies 
of incoming photons: (a) Q. = leV, and (b) Q, = 2eV. For 
each case, total spectral density g{Lij) and contributions due 
to the n~ -^n^ and (n ± 1)~ — >■ (n =p 1)+ modes are shown in 
the solid, dot-dashed, and dashed line, respectively. Intensity 
scale is the same on (a) and (b); values of the parameters 
used: v = lO^m/s, 71 = 0.4eV, B = lOT, and T = 0.012eV. 



portional to the inverse square of the incoming photon 
energy 57. Therefore, two peaks drawn with dot-dashed 
hnes are roughly four times smaller on the right figure. 
The spectral density of the (n ± 1)~ — (n =F 1)+ transi- 
tions is smaller by a further factor (^) in comparison 
to the n~ — > n"*" transitions. Hence, this contribution, 
shown with dashed lines, is close to zero on the right fig- 
ure, while on the left, the first of the two smaller peaks 
corresponding to symmetric transitions 2^ — > 0, — > 2+ 
and 3~ — 1, 1 — > 3+ is still visible in the total spectral 
density. Because of the contrasting polarization factors 
in Eq. (jlip . contributions of different modes, n~ — ^n"*" or 
(n ± 1 ) ~ — )• (n =F 1 ) , to the total spectral density could be 
separated using polarizers. If the polarizers were set as to 
collect only photons with circular polarization identical 
to that of the incoming photons, then the — >n+ con- 
tribution would be measured. However, if only the pho- 
tons with polarization opposite to the polarization of the 
incoming beam were detected, the (n =F 1)^ (n ± 1)^ 
contribution would be determined. 

Increasing the filling factor leads first to the 2^ 
and 3^ — > 1 transitions being blocked when LLs with n — 
and n = 1 are completely filled. Therefore, the height 
of the two corresponding (n±l)~— ^(?i=f1)+ peaks is 
halved (transitions — > 2+ and 1 — > 3+ are still allowed). 
Next to disappear are the first n~ — )■ peak, that is 
2^ — > 2+, and the remains of the first (n± 1)^ — (n^ 1) + 
peak, (due to the — 2+ transition) because of the filled 
LL 2+. Complete filling of each following Landau level 
results in the disappearance of the next — > 7i+ and 
(n ± 1)^ [n =F 1)+ peaks. 



IV. SUMMARY 

We presented a theory of inelastic scattering of pho- 
tons in bilayer graphene accompanied by the excitation 
of electron-hole pairs . Similar to monolayer graphene. 



the dominant scattering processes lead to the crossed po- 
larisation of in/out photons. Also, the selection rules 
in the presence of a magnetic field are found to be the 
same, with the n~ — >-n+ mode being the strongest. We 
estimate the intensity of one of the — > scattering 
processes to be /„-^„+ ^ 10^^'^ for Vl ^ leV photons 
in magnetic field B ~ lOT. The most recent theoretical 
prediction for monolayer graphene of the intensity of the 
phonon-induced G peakj^^ a well known Raman feature 
in carbon materialsf^l estimates Iq ~ 10""'^^. This result 
is only one order of magnitude greater than the intensity 
of a single rt^ — >■ peak. Therefore, spectral features 
of inter-Landau-level transitions in bilayer graphene pre- 
dicted in this paper may be observable experimentally. 

The electron Raman scattering in an external mag- 
netic field would complement infrared spectroscopy as 
it can give information about electronic excitations be- 
tween different pairs of Landau levels. The purpose of 
this paper was to identify and describe the dominant 
inter-Landau-level modes. However, additional correc- 
tions e.g., to the Raman-active magneto-exciton energies, 
due to the many-body effects neglected in the above con- 
siderations will be present in the spectra. Many-body 
corrections were, for example, observed in the infrared 
spectroscopy experiment performed in external magnetic 
fields Electronic Raman measurements could provide 
a test ground for some of the theoretical models proposed 
to account for these many-body effects. 
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Appendix: The contact interaction 

The contact interaction scattering amplitude (572,, il- 
lustrated using Feynman diagram in Fig. 1(d), results 
from the second term in the interaction Hamiltonian ([3]) 
and corresponds to the one-step Raman processes. It is 
characterised by operators d'^H/dpidpj and hence, the 
bilayer graphene Hamiltonian in ([1]) does not allow any 
contact interaction processes, as it includes only terms 
linear in the electron momentum p. The only contribu- 
tion to the contact interaction comes from the addition 
SH, which contains terms quadratic in p. However, such 

2 

contribution involves prefactor and therefore leads 
only to small corrections in the intensity of Raman scat- 
tering of photons with energy less than the band-width 
of graphene, ~ 670. In fact, the contribution to the scat- 
tering amplitude due to the contact interaction obtained 
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within the four-band model is 

2 fc2 2 

SlZ=^-^C-d; C={C^,Cy); (A.l) 

However, as we are interested only in the low-energy 
physics (small Raman shifts), only the top left block of 



the above matrix is relevant. This block contains an 
additional prefactor We see that the contact in- 

teraction can be neglected in comparison to the lead- 
ing terms in the effective scattering amplitude TZeS, as 

^ 71. ^ 1 

670 ^ 02 ^ n- 
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